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Harter,  H.  Leon  and  Moore,  Albert  H.,  local-maximum-likeliiiood  estima¬ 
tion  OF  THE  PARAMETERS  OF  THREE-FARAMETEn  LOGNORMAL  POPULA¬ 
TIONS  from  complete  and  CENSORED  samples,  Vol.  61,  No.  315  (Septem¬ 
ber,  10G5),  842-Sol. 

'I’ll ('80  correctioiiH  were  supplied  by  the  authors: 

I*.  842.  lloth  authors  arc  employed  at  Wright- Patterson  Air  Force  Base,  but 
the  second  author  is  with  tiic  Air  Force  Institute  of  Technology,  not  the  Aero¬ 
space  Research  Laboratories. 

P.  S45.  On  the  fourth  line  from  the  bottom,  the  lower  limit  of  the  integral 
should  be  4r+i. 

P.  840.  On  the  Inst  lino  of  equation  (3.1C),  the  sign  preceding  4r+i/(4r.»j)/5i 
should  be  positive. 

I*.  S47.  In  Table  1,  the  entries  in  the  second  column  from  the  right,  first  line, 
and  the  third  column  from  the  right,  third  and  fifth  lines,  should  be  —0.111959, 
-0.352750,  and  -0.513932,  respectively. 

P.  850.  In  Table  2,  the  last  two  column  headings  in  the  first  section  of  the 
lower  half  of  the  table  should  be  Cov  (fi,  t )  and  Cov  (&,  t),  respectively;  the 
value  of  Cov  (fl,  &)  when  N  =  200,  Q1  —0.00,  Q2«*0.0  should  be  —0.0010. 
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The  natural  logarithm  of  the  likclihniiil  function  in  written  down  for 
the  ni  -  r  order  statistics  remaining  after  censoring  the  n  —  m  largest  and 
the  r  smallest  observations  of  a  sample  of  sire  n(0<r<m<n)  from  a 
three-parameter  lognormal  population.  Its  first  partial  derivatives 
with  respect  to  the  parameters,  when  equated  to  zero,  yield  the  likeli¬ 
hood  equations,  and  the  negatives  of  its  second  partial  derivatives  with 
respect  to  the  parameters  are  the  elements  of  the  information  matrix. 
Algebraic  solution  of  the  likelihood  equations  is  impossible,  so  it  is 
necessary  to  resort  to  iteration  on  an  electronic  computer.  The  iterative 
procedure  proposed  is  applicable  to  special  cases  in  which  one  nr  two  of 
the  parameters  are  known  as  well  as  to  the  most  general  case  in  which 
all  three  parameters  are  unknown.  A  modification  of  the  procedure  al¬ 
lows  circumvention  of  a  certain  anomaly  which  sometimes  occurs  in 
maximum-likelihood  estimation  of  the  parameters  of  a  three-paramoter 
lognormal  population  from  small  samples.  The  information  matrix  is 
inverted  to  obtain  the  asymptotic  variances  and  covariances  of  the 
local-maximum-likelihood  estimators,  which  are  tabulated  for  various 
values  of  the  censoring  proportions  q i  «■  r/n  from  below  and  71  —  (n  —  m)/ 
n  from  above.  Ilcsulls  arc  reported  of  a  Monte  Cnrlo  study  conducted  to 
check  the  validity  of  the  asymptotic  variances  and  covariances  and  their 
applicability  to  samples  of  moderate  size. 

1.  INTKODUCTION 

In  TtiEtu  Ixmk  011  the  lognormal  distribution,  Aitchison  and  Brown  (1957, 
pp.  37-05)  have  devoted  two  entire  chapters  to  estimation  problems,  one 
each  for  the  two-parameter  and  three-parameter  distributions.  They  have 
given  a  comprehensive  summary  of  efforts  up  to  that  time  to  cstimato  the 
parameters  of  a  lognormal  population  by  the  method  of  maximum  likelihood, 
the  method  of  moments,  the  method  of  quantiles,  the  graphical  method,  and 
mixed  methods.  The  problem  of  estimating  the  parameters  of  a  two-parameter 
lognormal  population  with  known  lower  bound  r  is  equivalent  to  that  of  es¬ 
timating  the  parameters  of  a  normal  population,  which  has  been  considered  by 
a  number  of  authors.  Harter  and  Moore  (1900)  have  summarized  and  extended 
the  contributions  of  others  to  the  solution  of  that  problem,  with  particular 
emphasis  on  the  method  of  maximum  likelihood,  and  have  proposed  an  itera¬ 
tive  procedure  for  obtaining  maximum-likelihood  estimates  of  the  parameters 
from  complete,  singly  censored,  and  doubly  censored  samples. 

Maximum-likelihood  estimation  of  the  parameters  of  a  three-parameter 
lognormal  population,  for  complete  samples,  has  been  investigated  by  Wilson 
and  Worcester  (1945),  Cohen  (1951),  Aitchison  and  Brown  (1957,  pp.  55-6), 
and  Hill  (19C3).  The  latter  has  explored  some  unusual  features  of  the  likeli¬ 
hood  function  of  the  three-parameter  lognormal  population  which  had  appar¬ 
ently  gone  unnoticed  by  earlier  investigators.  In  particular,  he  has  shown  that 
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then:  exist  paths  along  which  the  likelihood  function  of  any  ordered  sample 
.'  i,  •  •  •  ,  J  .,  tends  to  as  (r,  n,  a )  approaches  (xi,  —  *,  +  «),  where  r  is  tho 
threshold  parameter  and  u  and  a  are  the  mean  and  the  standard  deviation  of 
the  parent  normal  population.  This  global  maximum  of  the  likelihood  function 
leads  to  the  ridiculous  maximum-likelihood  estimates  f=xi,  £=  — »,  and 
a  ~  -|-  X  regardless  of  the  sample.  On  the  other  hand,  solution  of  the  likelihood 
equations  leads,  in  most  eases,  to  local-maxim..  11-likelihood  estimates  which, 
while  not  true  maximum-likelihood  estimates  according  to  the  usual  definition, 
are  reasonable  estimates  and  appear  to  possess  most  of  the  desirable  properties 
usually  associated  with  maximum-likelihood  estimates.  Exceptions  may  occur 
in  the  case  of  small  samples,  for  which  the  likelihood  function  may  have  no 
dearly  defined  local  maximum. 

Apparently  nothing  has  been  published  on  the  problem  of  estimation  for 
truncated  or  censored  samples  from  a  three-parameter  lognormal  population, 
t  hough  Aitchisonnnd  Hrown  (11)57,  pp.  88-91)  have  discussed  the  two-parame¬ 
ter  case,  'flic  present  paper  will  lie  devoted  to  local-maximum-likelihood  es¬ 
timation,  for  the  three-parameter  case,  from  singly  anil  doubly  censored  ns 
well  as  complete  samples.  An  iterative  estimation  procedure  for  use  on  a  high¬ 
speed  computer  will  be  given.  This  procedure  will,  in  most  coses,  converge  to  a 
point  where  the  likelihood  function  has  a  local  maximum.  If  the  sample  (after 
censoring,  if  any)  is  small,  convergence  may  be  slow  or  the  iterative  procedure 
may  take  oft  along  the  path  to  infinity  mentioned  above.  Even  in  the  latter 
case,  reasonable  estimates  can  be  obtained  by  a  modification  of  the  iterative 
procedure  similar  to  that  used  by  Harter  and  Moore  (19f>r»)  for  the  three- 
parameter  Clamma  and  Weibull  populations. 


* 


*.  THE  MKEMHOOl)  EQUATIONS 

Consider  a  random  sample  of  size  n  from  a  three-parameter  lognormal  pop¬ 
ulation  with  parameters  m,  0,  and  r  (the  location  or  threshold  parameter). 
I.et  .Ym  1,  •  •  •  ,  .Yw  be  the  ordered  observations  remaining  after  the  n—m 
largest  and  the  r  smallest  observations  have  been  censored.  The  joint  prob¬ 
ability  density  function  of  these  order  statistics  is  given  by 


/(*rn.  ■  ■  ■  a,  r) 

__  n! _  TJ _ _1 _ (  y  [ln(z<  —  r)  -  m]*' 

“  (h  -  m)  !r!  <r VMx,  -  r)  \  tZU  ~  2** 


t,  rin(T,t.i  -  r)  -  Ml)  ' 

V  L  .  Jf  V 

L  9  J/ 

(2.1) 


The  natural  logarithm  of  the  likelihood  function  is  given  by 

1  "* 

L  =  In  [a!/ (it  —  m)!r!j - (m  —  r)  In  2ir  —  (m  —  r)  liia  —  ln(z<  —  r) 

2  l—r+ 1 


1  "  , 

4-  E«!  +  (» 


m)ln  [l  -  Fir.)]  +  rlnF(*r+i). 


(2.2) 
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u  =  [luUi  -  r)  sj/<r,  F(z()  -  J  /(*)«'<,  and  /(*,)  =  (2ir)  1  cxp(-z*/2). 


The  likelihood  equation  are 


dl.  ! 

a 

OL  l 


,  ,  4  /(*»)  /(*,♦.) 

+  («  —  m)  -  —  r - 

l  -  F(zm)  FM 


da  a 


{  t  '■ 

I  l-r+l 


(  ,  ,  ,  V'  ,  ,  ,  ,  Zm]{Zm)  J,+  |/(rr+l 

l  ‘  1  -  FM  F(zf+1) 


/Ur+l)  1 

r - >  -  0, 

FM) 

Zmf{Zm)  Zr+lf(Zr+\) 


}-o. 


E  (■*■« —  r)  1  + 


Li  £  Jl- 

9  \  i«r-fl  J* i  T 


+  («  -  m) 


(/„-r)[l  -F(jm)] 


r _ /(*rf.) _ ) 

r  (.r,+  i  -  r)F(zr+,)j 


0.  (2.5) 


If  m  *»  re  and  r»(),  equations  (2.11)  ami  (2.4)  can  be  solved  for  n  and  <r  as  explicit 
functions  of  r,  yielding 

P  -  E  M*.  -  t)]/r,  (2.0) 


*  -  |/  E  IW*.  “  r)  ~  P)’/n. 


The  equation  (2.5)  cannot  he  solved  explicitly  even  if  m  — n  and  r>»0.  If 
censoring  occurs,  none  of  the  likelihood  equations  has  an  explicit  solution,  and 
it  is  necessary  to  resort  to  iterative  solutions.  The  details  of  an  interative 
procedure  for  solving  the  likelihood  equations  will  be  given  in  section  4. 

3.  ASYMPTOTIC  VARIANCES  AND  COVARIANCES  OP  ESTIMATORS 

The  second  partial  derivatives  of  L  with  respect  to  the  parameters  arc 

m  i  i . /(*-)  r  /(*->  i 

5?  _  _|_(m  -  r)  +  (n  -  m)  ~  -^-J 


JM 


*r+i  + 


/(«>♦.) ; 

F(z,+t). 


d'L  1 
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Estimation  up  tiihee-pahametek  loonohmai,  poi-ulationh 
/(*») 
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(3.3) 
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<h  =  >’/«.  ’It  =*  («  -  m)/n,  uml  p  =  I  -  9,  -  7,  =  (m  -  r)/n. 
As  n— * *>  (</,  and  q,  fixed),  ?r+,-*fr+,  where 

J  f(t)dl  -  }1,  *»,  — •  i„  where  f  f(l)dt  «  q,, 
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,_,+i 
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t':  £  !*./(-»•.  -  r)(  -»  f  e-<-+**>/l/(i)rf< 

.W+' 


+  <r)  -  F(i,+ 1  +  a)] 
4-  /(i«  +  <r)  -  /(2r+I  +  *)}, 


K  |)  |(<r;,  -  I)/(x«  -  r)1  j  -»  f  ,  -*<n"U0t  _ 


t  )/«)<« 


=  c5,#,'<‘,{(2ff*  +  Dim, +  2,) 

-  *’(*.„  +  2<r)|  +  +  2<r)  -  /( f,+>  +  2<r)]} . 

l  lic  elements  of  the  information  matrix  (multiplied  by  <r*/«)  may  be  written  as 

Jim  (-<,yn)£(av,/<V)  -  />  -  f(i»)/q2] 

+  /(4,i)|U,  +/(«m)/ji|  -  r»,  (3,1I) 

lim  (-»*/«) /••<**/,/,>„*>  -  2 p  -  i4(im)  -  -  /(ft.)/,,,] 

(3.12) 

4-  2r+l/(2r+l)  +  2r+l/(2r+l)  [ir+ 1  + /(2 r+l)/(7l]  =>  I'**, 

Jim  (-ffVn)A’(d3/y«lr*)  -  r*'*’-*' { („*  +  l)[f’(2„  +  2<r)  -  F(2,+I  +  2a)\ 

+  »!/(«-  +  2<r)  -  /(«  r+1  +  **)]} 

-  [«-  -  fU~)/q,  +  <r]  (  I3) 

+  e“,,“**'*<«,/(fr+l)(fr+,  +  /(ft+1)/7,  +  a]  «  e”, 

lim  (-e'/n)fC(d'-L/<W<,)  ~  -  /(ft.)  -  V(*-)[ft.  -  /(U/f»|  +  /(ft+I) 

*"**  (3.14) 

+  £rH/(f,+  |)[ir+l  +  /(ir+l)/l?l]  “  I’”, 

Jim  i-V/nWd'L/dtfr)  «  +  <r)  -  F(ft+,  +  <r)] 

+  »■-'»+«♦, •>/(«r+I)f4r+,  +  /(ft+,)/<7tl  -  e". 
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(3.15) 


+  /(&.  +  ")  -/(*,+ 1  +  e)| 

-  «-<*♦•-'’/(«*)  [£  -  l^«.)/?,  -  1] 


(3.18) 


+  e-<»+**'*,»»/(ir+,)(4i  -  -  lj  -  #**. 

The  asymptotic  variance-eovariancc  matrix  for  the  estimators  jft,  *,  and  f  is 
then  (<r*/n)[p(yl,  where  fe(/J«=  [»*']-'.  The  reader  will  note  that  t>",  v*,  and  e>‘ 
are  independent  of  both  a  and  <t,  while  b*\  e”,  and  e**  are  dependent  upon  both. 
The  dependence  upon  a  is  only  through  the  factor  <r*  in  r1*  and  i/**  and  the 
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factor  c  •*  in  t,M,  while  the  dependence  upon  a  is  more  complicatetl.  It  follows 
that,  among  the  elements  of  the  inverse  matrix  [t*,/ 1  =*  f'n,  pt»,  and  vu 

are  independent  of  g,  in,  ami  e.»  depend  upon  a  only  through  the  factor  ef,  and 
i  is  depends  upon  g  only  through  the  factor  c1*. 

The  computation  of  the  elements  v 11  of  the  information  mfttri.v  ^multiplied 
hy  o-/n)  and  the  inversion  of  this  matrix  and  its  aubmatriccs  and  multiplica¬ 
tion  hy  o1  to  obtain  the  coelhcicnts  of  \/n  in  the  asymptotic  variances  and 
covariances  were  performed  on  the  IBM  1020  computer.  The  resulting  co- 
eHicients  of  |/«  in  Var  (/J),  Vnr  (&),  Var  (#),  Cov  (ft,  &),  Cov  (ft,  f),  and  (’ov 
(d,  t)  are  given  in  Table  I  for  g=l,  <r  =  2  with  censoring  proportions  71  =  0.00, 
0.01 , 0.02  and  r/i  =■-  0.0,  0.5,  together  with  the  coeflicients  of  1  /n  in  the  variances 
and  covariances  when  one  parameter  is  known  and  in  the  variances  when  two 
parameters  are  known.  For  <r=^2  and  other  values  of  g,  the  coeflicients  can  be 
obtained  from  the  tabular  values  in  Table  I  by  multiplying  by  the  proper 
exponential  function  of  g,  but  for  other  values  of  n  the  computations  must 
be  made  afresh.  The  coefficients  tabulated  are  those  of  the  local-maximum- 

TABLE  t.  COEFFICIENTS  IN  ASYMPTOTIC  VARIANCES  ANI)  COVA1UANCK8 
OI’  MAXIMUM-LIKELIHOOD  KSTIMATOltS  OF  PARAMETERS  g,  *,  AND  r  OF 
LOC, NORMAL  POPULATION  WITH  #-2  FROM  SAMPLE8  OF  SIZE  .V  WITH 
PROPORTIONS  qi  CENSORED  FROM  BELOW  AND  Q2  FROM  ABOVE 


Ul 

^2 

JV  VAR  <£> 

| 

.V  VAR  (S) 

N  VAR  It) 

.VCOV  d,0) 

at  rov  (£.?) 

at  COV  <#.t) 

exp  <j*-m 

EXP  (»— 4> 

EXP  1.-4) 

0.00 

0  0 

4  015152 

2.000005 

0.927274 

-0.030904 

-10.111000 

0.223919 

0.00 

0.5 

0.070045 

5.254124 

0.950419 
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likeliliood  estimates  resulting  from  the  iterative  estimation  procedure  which 
will  be  discussed  in  Section  4,  not  those  of  the  global  maxima.  This  fact  is 
verified  hy  the  results  of  a  Monte  Carlo  study,  which  will  bo  reported  in  Sec¬ 
tion  5,  comparing  variances  and  covariances  of  estimates  obtained  from  500 
random  samples  of  each  of  various  sires  and  various  degrees  of  censoring  with 
those  given  hy  the  asymptotic  formulas. 

4.  ITERATIVE  ESTIMATION  PROCEDURE 

The  procedure  for  iterative  estimation  on  a  high-speed  computer  involves 
estimating  the  three  parameters,  one  at  a  time,  in  the  cyclic  order  jt,  a,  r. 
omitting  any  assumed  to  l»c  known.  Assuming  that  the  first  m  order  statistics 
of  a  sample  of  size  n(m<n)  are  known,  except  for  the  first  rn(0$r#<m  —  ;>, 
where  p  —  I,  2,  or  3  is  the  number  of  parameters  to  be  estimated),  one  starts 
by  setting  r « r0.  One  then  chooses  initial  estimates  for  the  unknown  parame¬ 
ters. 

At  each  step,  the  rule  of  false  position  (iterative  linear  interpolation)  is  used 
to  determine  the  value  (if  any)  of  the  parameter  then  being  estimated  which 
satisfies  the  appropriate  likelihood  equation,  in  which  the  latest  estimates  (or 
known  values)  of  the  other  two  parameters  have  been  substituted.  For  9 
<xr,u,  one  ran  always  find  estimates  it  (finite)  and  a  (finite  and  positive)  in 
this  way.  In  estimating  r,  however,  one  may  find  that  no  value  of  r  in  the 
permissible  interval  satisfies  the  likelihood  equation  (2.5).  In  such 

cases,  the  likelihood  function  is  monotone  increasing,  so  that  f  =  xr,+i.  As 
9— p-4—  «  and  #—*+*,  so  that  the  estimation  is  proceeding  along  a 
path  to  the  globul  maximum  [see  Ilill  (19(53)].  When  this  occurs,  which  is  not 
unusual  when  the  available  sample  (after  censoring,  if  any)  is  small,  it  is  still 
possible  to  obtain  reasonable  estimates  by  a  slight  modification  of  the  proce¬ 
dure.  The  modification  entails  censoring  x,,+i,  the  smallest  observation  not 
previously  censored,  and  any  others  equal  to  it,  thus  increasing  r  from  r#  to 
cu+ri,  when>  ri>l  is  the  number  of  observations  censored  at  this  point.  Sub¬ 
sequently,  j„»i  plays  no  role  in  the  estimation  procedure  except  as  an  upper 
bound  on  9.  Now  the  likelihood  function  is  bounded  and  finite  estimates  (l  and 
b  are  obtained.  Iteration  continues  until  the  results  of  successive  steps  agree  to 
within  assigned  tolerances  (say  10“*)  or  for  a  specified  number  of  steps  (say 
550). 

Use  of  the  modified  procedure  calls  for  distributions  of  estimates  which  are 
conditional  on  the  necessity  for  censoring  the  smallest  previously  uncensored 
observation (s).  backing  these,  one  may  assert  that  the  asymptotic  variances 
of  the  estimators  for  the  “sometimes  censor”  procedure  are  lioimded  below  by 
those  for  the  “never  censor"  procedure  and  alxvve  by  those  for  the  “always 
censor”  procedure. 


S.  MONTE  CARLO  STUDY  OF  ESTIMATES 

In  order  to  check  the  validity  of  the  asymptotic  variances  and  covariances 
determined  in  Section  3  and  their  applicability  to  samples  of  moderate  size,  a 
Monte  Carlo  study  was  carried  out  on  the  IBM  7094  computer.  Five  hundred 
pseudo-random  samples  each  of  sizes  50,  100,  and  200  from  a  lognormal  pop- 
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iilutiun  with  parameters  it  — 4,  ’<r  —  '2,  and  r-  10  were  generated  in  the  com¬ 
puter.  The  iterative  procedure  described  in  Section  4  was  used  to  estimate  all 
three  parameters,  also  every  pair  of  parameters  and  every  single  parameter, 
the  known  values  being  substituted  for  the  parameters  not  l>eing  estimated.  In 
t  lie  ease  of  samples  of  size  100,  this  was  done  not  only  for  the  complete  samples 
hut  also  for  the  samples  with  proportions  f/,  -0.01  censored  from  below  and 
</•_'  0.5  from  above,  both  singly  and  in  combination.  The  unmodified  proce¬ 
dure  converged  to  within  the  assigned  tolerances  of  10**  in  the  specified  550 
steps  or  fewer  without  exception  for  the  uncensored  and  the  singly  censored 
samples.  Among  the  doubly  censored  samples,  however,  there  were  40  in 
which  it  was  necessary  to  resort  to  the  modified  proved um  (censoring  the  second 
order  statistic)  in  order  to  estimate  all  three  parameters  simultaneously  and  23 
in  which  it  was  necessary  to  do  so  in  order  to  estimate  a  and  r  simultaneously 
with  n  known.  More  severe  censoring  or  the  use  of  smaller  samples  can  lead 
to  slower  convergence  of  the  iterative  procedure  as  well  as  to  the  need  to  resort 
to  the  modified  procedure  for  a  greater  proportion  of  cases. 

The  means,  variances,  and  covariances  of  the  estimates,  based  on  500  Ham- 
pics,  are  given  in  Table  2.  For  the  doubly  censored  case,  Table  2  gives  two  sets 
of  results,  one  excluding  and  the  other  including  the  coses  in  which  it  was  neces¬ 
sary  to  resort  to  the  modified  procedure.  A  comparison  of  the  means  in  Table  2 
wit  h  the  population  parameters  »  -  4,  a  -  2,  and  r  =  10  and  a  comparison  of  the 
variances  and  covurinnces  in  Table  2  with  the  asymptotic  variances  and  co- 
variances  obtained  by  dividing  the  coefficients  in  Table  I  by  the  sample  size 
leads  one  to  the  following  tentative  conclusions: 

(1 )  The  estimator  A  from  complete  samples  has  a  negative  bias  which  is  ap¬ 
proximately  proportionnl  to  ir\  the  reciprocal  of  the  sample  gizc.  The  bias  of 
A  appears  to  be  unaffected  by  knowledge  of  <r,  but  it  is  small  or  non-existent 
if  r  if  known.  The  bioH  remains  negative  and  increases  in  absolute  value  for 
moderate  censoring,  but  may  be  positive  in  cases  of  severe  censoring — see  con¬ 
clusion  (4). 

(2)  The  estimator  t  has  a  positive  bias  when  r  is  unknown  and  a  negative 
bias  when  r  is  known.  The  magnitude  of  the  bias  is  roughly  proportional  to 
m"\  and  appears  to  tie  unaffected  by  knowledge  of  it.  It  is  increased  by  censor¬ 
ing,  especially  from  Mow. 

(3)  The  estimator  9  has  a  positive  bias  which  is  closely  proportional  to 
n~K  The  magnitude  of  the  bias  appears  to  l>e  unaffected  by  knowledge  of  it 
and/or  a  and  by  censoring  from  above,  but  if  is  markedly  increased  by  cen¬ 
soring  even  a  single  observation  from  below. 

(4)  Cases  in  which  it  is  frequently  necessary  to  resort  to  the  modified  pro¬ 
cedure  arc  characterized  by  fairly  large  positive  biases  of  all  the  estimators, 
especially  &  and  9.  Tho  biases  arc  only  moderately  increased,  however,  by  in¬ 
cluding  rather  than  excluding  the  instances  in  which  resort  to  the  modified 
procedure  is  necessary. 

(5)  The  variance  of  9  for  samples  of  moderate  size  is  much  larger  than  the 
value  given  by  the  asymptotic  formula  (sec  Table  1).  The  excess  over  the 
asymptotic  value,  however,  appears  to  be  proportional  to  n-*,  whereas  tho 
asymptotic  value  itself  is  proportional  to  n-1.  Thus,  for  sufficiently  large  sam- 
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TAD  LB  'i.  MKANS,  VARIANCES,  AN  I'  COVARIANCES  OK  KSTI  MATES  OF 
PA  II A  METE  IIS  FROM  500  RANDOM  SAMPLES,  EACH  OF  SIZE  N,  WITH 
PROPORTIONS  Ql  CENSORED  FROM  1IELOW  AND  Q2  FROM  AROVE,  DRAWN 
FROM  LOGNORMAL  POPULATION  WITH  PARAMETERS 
m“”I,  fm‘£,  AND  f  «  10 
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*  Keeiilt*  giveti  ocj  thi*  liw  tr t  for  »ll  500  cuw,  including  those,  omitted  on  the  line  •hove,  which  required  use 
of  the  modified  procedure  (censoring  tbe  eecood  order  •tetietic),  eo  the!  Ql  *0.02. 
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plus,  tho  excess  would  become  negligible  in  comparison  with  the  asymptotic 
value.  The  excess  appears  to  be  decreased  little  if  any  by  knowledge  of  a  and 
only  slightly  by  knowledge  of  <r.  'The  excess  appears  to  l>c  affected  very  little 
by  censoring,  but  for  censoring  of  even  a  single  observation  from  below  the 
asymptotic  value  increases  markedly,  so  that  the  excess  becomes  smaller  by 
comparison. 

(<>)  When  r  is  known,  the  variances  of  p  and  &  are  in  close  agreement  with 
the  values  given  by  the  asymptotic  formulas.  When  t  is  unknown,  the  vari¬ 
ances  of  (i  and  it  and  their  covariances  with  t  are  somewhat  larger  than  the 
values  given  hy  the  asymptotic  formulas.  The  excess  over  the  asymptotic 
value  is  greater  for  a  than  for  p,  as  one  would  expect  from  the  fact  that  6  is 
more  strongly  correlated  with  f  than  is  p. 

(7)  For  cases  in  which  it  is  frequently  necessary  to  resort  to  the  modified 
procedure,  tho  variances  and  covariances  including  and  excluding  instances  in 
which  the  modified  procedure  is  actually  used  do  not  appear  to  differ  sys¬ 
tematically  from  each  other  or  from  those  given  by  the  asymptotic  formulas, 
except  that  Cov  (&,  f)  is  unexpectedly  small. 

The  authors  wish  to  thank  the  Editor  and  the  referees  for  a  number  of  help¬ 
ful  suggestions. 
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